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The constituent quark number scaling of elliptic flow is studied in a non-equilibrium hadronization 
and freeze-out model with rapid dynamical transition from ideal, deconfined and chirally symmetric 
Quark Gluon Plasma, to final non-interacting hadrons. In this transition a Bag model of constituent 
quarks is considered, where the quarks gain constituent quark mass while the background Bag-field 
breaks up and vanishes. The constituent quarks then recombine into simplified hadron states, while 
chemical, thermal and flow equilibrium break down one after the other. In this scenario the resulting 
temperatures and flow velocities of baryons and mesons are different. Using a simplified few source 
model of the elliptic flow, we are able to reproduce the constituent quark number scaling, with 
assumptions on the details of the non-equilibrium processes. 
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I. INTRODUCTION 

It has been observed that the momentum distribu- 
tion of particles created in heavy ion collisions is not 
azimuthally symmetric in the plane perpendicular to the 
beam direction. This dominant asymmetry, the elliptic 
flow, is the result of several factors, the main one being 
the anisotropy of the initial configuration after the colli- 
sion due to the non-zero impact parameter. The elliptic 
flow is typically characterized using the second coefficient 
of the Fourier expansion of the momentum distribution, 
V2- It was found in experiments [IJ that the V2 parame- 
ter as a function of the transverse momentum, pt, scales 
with the number of constituent quarks, n^q, in a hadron. 
That is, if the V2{pt) curves are re-scaled according to 
the constituent quark number of the considered hadron 
species, and V2/ncq is plotted as function of pt/ricq for 
each type of hadron with mass , the curves will coin- 
cide. Later results showed that the scaling is more precise 
if V2 is plotted as a function of the transverse kinetic en- 
ergy, Et — ^ m\ + pI — rrih, instead of the transverse 
momentum pf 

This experimentally found scaling law is remarkably 
simple, and indicates that the elliptic flow develops be- 
fore the quarks recombine into hadrons. Therefore, un- 
derstanding the factors influencing the measurable V2 can 
provide information about the state of the Quak Gluon 
Plasma (QGP). 

The fluid dynamical (FD) model describes the dynam- 
ical development of the QGP from the (already thermal- 
ized) initial state until the break-down of the equilibrium, 
where first the chemical equilibrium among quarks and 
anti-quarks ceases. Initially the plasma has only two fla- 
vors, u and d, and thus the FD model assumes two fla- 
vors. However, the flavor equilibration in the plasma is a 
rapid process with a timescale of the order of 1 fm/c, so 
by the end of the FD development a flavor equilibrium 



among three flavors, u, d, and s, is reached. Our usual 
FD evolution does not take into account this chemical 
change of flavors with an additional rate equation, in- 
stead this process is taken into account at the flnal stage 
of the FD development, when the subsequent EoS is al- 
ready assumed to have three equilibrated flavors. En- 
ergy and momentum conservation and the requirement 
of non-decreasing entropy is enforced in the transition, 
from the ideal QGP state to the state where quark and 
anti-quark numbers are frozen out. As a consequence, 
the mass change of the quarks starts in the initial QGP 
with two flavors, and we use this approximation to esti- 
mate in Section II. B the final boundary where the initial 
FD stage of the evolution ends, with light quark masses. 

We consider a model of hadronization and investigate 
the constituent quark number scaling (QNS) of the V2 pa- 
rameter. In this model, a gas of quarks and anti-quarks 
expands in a background field represented by the Bag 
constant B, which depends on density and temperature 
of the expanding fireball. Initially this S-field includes 
the energy of the deconfined perturbative vacuum and of 
the gluon fields. As the system expands the deconfine- 
ment starts and the average B decreases. Furthermore, 
as the chiral symmetry breaking starts, the quarks gain 
mass. The quark mass is calculated as a function of the 
temperature and density of the matter. The quark gas 
expands rapidly while the quark mass increases. This 
process can be considered as a simple representation of 
the breaking chiral symmetry and deconfinement in a dy- 
namical transition crossing the Quarkyonic phase [2]. 

The point during the expansion when the quarks re- 
combine into hadrons is determined from the condi- 
tion that the average hadron energy is equivalent to 
(1.0 — 1.1) GeV, as found from the systematics of experi- 
mental data J3j. At recombination, the thermal and flow 
equilibrium between particles is broken. 

Finally, the V2 parameter is determined using simple 
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two- and three-source models of the cUiptic flow and the 
particle distributions obtained from the hadronization 
model. 

The paper is organized as follows: In Section [IT] 
the density and temperature dependence of constituent 
quark mass is considered. The energy, entropy, momen- 
tum and chemical potential of each individual source and 
of the total source is given in Section [Till In Section |TV] 
we describe the applied non-equilibrium model of rapid 
hadronization. Details of calculation of elliptic flow pa- 
rameter V2 are presented in Section [V] A summary is 
given in Section [VT} 

II. CONSTITUENT QUARK MASS 

To present our arguments we consider the Nambu- 
Jona-Lasinio model (NJL) which is motivated by 
Quantum Chromodynamics (QCD) and is basically a 
purely quark-quark-interaction theory; a comprehensive 
overview about NJL has been presented in [4]. The ex- 
tended NJL-model contains three types of quark-quark- 
interaction and is given by the Lagrangian 

f—u,d,s 
f—u,d,s a—0 

E E[fe7.A"'//)%fe75 7.A'^'Z/)'; 

f—u.d,s a—1 

+ ^ [det {qjil + -i^)qf) + det - 75)9/)] , (1) 

where — q^jla, 7'' are the Dirac matrices, m/ is the 
current quark mass of flavor /, Gs is the coupling con- 
stant of the scalar-current interaction, Gy is the cou- 
pling constant of the vector-current interaction terms 
and Gd is the coupling constant of the determinant al 
flavor-mixing term (determinant in flavor space). The 
Gell-Ma nn m atrices are A° where a = 0, 1, 2, . . . , 8 with 
A° = ^2/3 /, where / is the unit matrix. This La- 
grangian is used to derive the relation between con- 
stituent quark masses and chiral condensates, both for 
vacuum and for finite densities and temperatures. 

A. Constituent quark mass in vacuum 

First, we consider the constituent quark mass in vac- 
uum. For three flavors, SUf{3), the constituent quark 
masses in vacuum = Mf {ub = 0, T = 0) where / — 
M, d, s, are related to the ciral condensates in vacuum as 
follows [SH7] 

M° = - 2 Gs {uu)o - Go (dd)o {ss)o , (2) 



MO = ma-2Gs {dd)o - Go {uu)o {ss)o , (3) 
M° = - 2 Gs {ss)q - Gd {uu)o (dd)o . (4) 

Here u — jo and so on, m„ = 5 MeV, nid — 9 MeV and 
nig = 130 MeV are the current quark masses of the u- 
quark, d-quark and s-quark, respectively. Typical values 
of the chiral condensates in vacuum are given by, e.g. 
[SI [To] and references therein: 

{qq)o = (?Iu)o = (drf)o = -(0.225 GeV)3, (5) 
(ss)o = 0.7 {uu)o . (6) 

Typical values for the coupling constants are G5 = (15 — 
20) GeV-^ Gv ^ 0.5 Gs and Gd = -(160-240) GeV"^ 
[llin|[7]. We note that in relations ([2]) - Q the magnitude 
of terms proportional to Gd are small compared to the 
terms proportional to G5. 

B. Constituent quark mass in a hot and dense 
medium 

Let us now consider the case of quarks in a hot and 
dense medium, at early times of FD evolution. At this 
stage we assume to have two flavors. Then, the relations 
([2]) - Q for two flavors SUf{2) are reduced to 

Mf = mf-2Gs {qq)ns,T. (7) 

The suffix nB,T at chiral condensate denotes the Gibbs 
average over hadron states and meson states [9l [TTI - H3] 
and references therein. Note that for isospin-symmetric 
matter there is no difference in the density and temper- 
ature dependence of u-quark and d-quark condensate, 
that means {qq)nB,T = {uu)nB,T ~ {dd)nB,T- In the 
limit of high densities and temperatures the constituent 
quark mass, Mf, approaches the current quark mass, m/. 
The well-known model-independent linear density depen- 
dence of the chiral condensate has widely been applied in 
many investigations, e.g. [7| [51 [TH [H] . The temperature 
dependence of the chiral condensate has been determined 
up to order 0{T^) in [T6^. In order to determine the den- 
sity and temperature dependence of the chiral condensate 
we follow the arguments of Refs. [71 [TH [Tl], where the 
first leading terms in the low-density low-temperature 
expansion have been obtained: 

{qq)nB,T = {qq)o 

( 3(7 A \ 

l,^ ^/f^''^" 871 " 38471 " 28871 

(8) 

The temperature and density dependence of chiral con- 
densate has been plotted in Ref. [14]. The baryonic den- 
sity in terms of quark degrees of freedom is given by 

'^B = ^ E ("/ ~^/) ' (9) 

f=u,d 
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where n/ (uf) is the quark (anti-quark) density. The 
baryonic density in a given volume V is related to the 
conserved baryon number Nb by = note the 

relation he = 197.3 MeVfm. For the logarithmic scale 
we take ~ 300 MeV, for details see [T6j. The pion 
mass in vacuum is m^r = 138 MeV and the pion decay 
constant in vacuum is = 93 MeV. The numerical 
value of the quark-sigma-term is aq — 15 MeV (see e.g. 
Ref. [7]), which is three times smaller than the nucleon- 
sigma-term ajy = A5 MeV. 



Mf [MeV] 




n [fm-"] 



FIG. 1; Temperature and density dependence of the con- 



stituent quark mass Af/ according to Eq. ( 10 1 



By combining Eq. ^ with Eq. ([s]) , we obtain the expres- 
sion for the in-medium mass of constituent-quarks: 



Mf =mf - 2Gs{qq)n 



3gq 



8fl 



rp4 

384/4 



288/6 T 



(10) 

The temperature and density dependence of the con- 
stituent quark mass is plotted in FIG. [T] According to 
Eq. (10 1, for sufficiently high temperatures and densities 



the constituent quark mass will coincide with the current- 



quark-mass mf. Of course, the applicability of Eq. (10) 
is restricted to low densities and temperatures. In any 
case, the most upper limit for density and temperature is 
given by the condition that constituent quark mass has 
to be positive. 



III. EXPANDING QUARK-ANTI-QUARK BAG 

The Bag consists of quarks, anti-quarks, and the gluon 
fields. Non-perturbative effects in the QCD vacuum and 
the energy density of gluon fields are described by the Bag 
constant, B. In our model, we subdivide the bag into N 
individual cells, each of which moves with an individual 
fiow velocity Vi. 



A. Total energy of the bag 

The total energy of the bag in the center-of-mass frame 
of the colliding nuclei is given by the volume integral of 
the " 00" component of the energy momentum tensor of 
each cell i, T/"" = {eJ + P') uf^u";; - gf"" . Considering 
that our EoS is given as a sum of the energy and pressure 
of the ideal quark and anti-quark gas plus the uniform 
bag energy density B, that means e* = e* + e* -f- _B and 
pi ^ pi _,_ 75* _ tiig total energy of aU cells of the 
bag in the center-of-mass frame of the colliding nuclei is 
given by 



N 



total 



(7f (e 



FM - 



N 

E 



V'P\ (11) 



where the sum runs over the number of all cells of the 
bag. Here, we assume a uniform bag-field energy den- 
sity B over the whole volume, V; B depends on density 
and temperature of the bag, thus on the time of the evo- 
lution of the fireball. We have used the notation e* for 
the invariant scalar, rest-frame energy density, and P* 
for the scalar pressure, V'' is the volume of cell i in the 
center-of-mass frame, and 7^ = 1 / -^Z 1 — v| where is 
the 3-velocity of cell i. In order to determine the rest- 
frame energy or proper energy of the cells, we assume 
a Boltzmann-Jiittner distribution [T7] for the particles 
inside each individual cell: 



1 



{2 7rhy 



exp 



T, 



(12) 



and / (Px J (Pp fi{x,p) = Ni is the normalization of the 
Jiittner distribution, where is the number of particles 
of a given type inside cell i under consideration. The 
proper energy density is a Lorentz invariant quantity by 
definition and, therefore, can be evaluated in any frame, 
e.g. in the local rest frame. Accordingly, we obtain the 
following expression for the energy density of quarks and 
anti-quarks for each cell of the bag: 



el + el + B.^ 



^ E Af|T,exp( 

f=u,d \ 



T, 



3K. 



^ E M]T,^.V 



S=uA 

Mf 



T, 




(13) 



(14) 



(15) 



Here, the density and temperature dependent constituent 



quark mass Mf is given in (10 1, and and are the 



baryonic density and temperature, respectively, of the 
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cell i. The initial value of the Bag constant is abbre- 
viated by Bq, and we use the numerical value Bq — 
(198 MeV)^ — 0.2 GeV fm~'^; see also considerations in 
Ref. [ig. 



B. Total entropy of the bag 

The total entropy in the center-of-mass Lorentz frame 
(CM) is given by 



N 

5total = ' 



(16) 



where the sum runs over all cells of the bag. The invariant 
scalar entropy density, s*, can be evaluated in any frame. 
The total entropy density consists of entropy density of 
the quarks and anti-quarks, and for the assumed Jiittner 
distribution is given by: 



(17) 



2^2 E ^^/-P( 

f=u,d \ 



MfKi 



Mf Ki 



f=u,d \ ' 




(18) 



(19) 



In what follows, we will take the adiabatic limit, which 
implies the total entropy to be a constant. We note that 
there is no Bag constant in the expression for the total 
entropy because the Bag-field is uniform, i.e. the entropy 
of Bag-field vanishes. 



C. Total momentum of the bag 

The total momentum of the bag in the center-of-mass 
Lorentz-frame is given by 

N 

O^Ptotai - E i^f {e' + P')v., (20) 



where the sum runs over all cells of the bag. The scalar 
pressure of each individual cell of the bag is given by 



= Pl + P^-B, 



(21) 



^-iE ^^exp(l) 

f=u,d \ J 



P. = 



1 

2^ 



E MfT^e^p 



f=u,d 




(22) 



(23) 



We note, that there is no Bag constant in the expression 
for the total momentum ( 20 ) because it ca ncel s out in 
the sum e' -I- according to Eqs. (13 1 and (21). 



D. Chemical potential of quarks and anti-quarks 

N) 

Assuming a Jiittner distribution, the density = 
of quarks having flavor / and the density of anti-quarks 

_ n\ . 

ny — having flavor / in a given cell i is determined 
by the following equations: 



Mf T, exp 



1 

2^ 
-^M^T.exp 




(24) 
(25) 



In mechanical (pressure), thermal and chemical equilib- 
rium, the number of quarks depends only on baryonic 
chemical potential fig and temperature T. Especially, if 
the quarks and anti-quarks are in chemical equilibrium 
with each other, then fiq — ^b/^ and Ji^ = —j-iB/'^- In 
order to determine the baryonic chemical potential we 
have to use the total baryonic number density ub given 
in Eq. ([9|. Then, according to the expressions (24) and 
(25), we obtain a definition of the chemical potential /x^ 
oFcell i as 



1 

3^ 



Ti sinh 



3T, 



E 

f=u,d 



Ml 
T, 



,(26) 



where we have used sinh a; — ^ (e^ — e^^). We assume 
that these thermodynamic relations hold in the ideal 
QGP until the final stage of the expansion starts. Then 
the net baryon number in each cell N^^ and the num- 
ber of quarks and anti-quarks, A^* ~ (n^ + n^j)Vi and 

— respectively, are given. While Vi and 

thus and rfj are changing during the expansion, the 

numbers of quarks and anti-quarks, and A'^* , are fixed 
at this point in each cell and do not change afterwards. 
This is the chemical freeze out in the model. 

The volumes, baryonic and quark densities and tem- 
peratures of each individual cell change during the ex- 
pansion. The chemical equilibrium between the baryon 
density and quark- and anti-quark densities breaks down 
at this point. Consequently, the chemical potentials for 
the single quarks are related to the single quark (and anti- 
quark) densities only and are not connected to the net 
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baryon density. Although we need the chemical poten- 
tials for the calculation of the energy density, pressure 
and entropy, these must be calculated from the given 
quark an d anti -qu ark densities in each cell, by using 
Eqs. (24 1 and (251, which determine and /I^. The 
net baryon nuniber is conserved and is determined 
by Eq. ^ using Eqs. ( |24[) a nd ( |25[ ). Only in case of 
chemical equilibrium Eq. ( |26[ ) is used to determine the 
chemical potential /i^. That is, the Ns/Nq and Ns/Nq 
ratio freezes out first at the initialization of our model cal- 
culation for the final expansion stage, where quarks start 
to gain mass, the background field B starts to disappear, 
but the constituent quark numbers do not change. The 
flow does not evolve, the cells expand while coasting with 
negligible pressure. Still considerable changes happen re- 
garding the masses and temperatures of the system. 



IV. NON-EQUILIBRIUM EXPANSION STAGE 

A. Rapid hadronization hypothesis 

A first order phase transition in chemical and thermal 
equilibrium with homogeneous nucleation would take a 
long time 19 , longer than two particle correlation mea- 
surements indicate. If the homogeneous nucleation can- 
not support the required rapid transition then the tran- 
sition becomes delayed and freeze out and hadronization 
will happen rapidly and simultaneously from a super- 
cooled QGP. Thus, a rapid process must be out of equi- 
librium, at least of chemical equilibrium [20]. Possible 
detailed mechanisms of this out of equilibrium transition 
are addressed recently in several works [211 ES] . In the 
framework of the present model this transition is repre- 
sented by the mass gain and coalescence of constituent 
quarks, and the simultaneous disappearance of the Bag- 
field. These two processes are treated phenomenologi- 
cally, while we enforce conservation laws and our model 
constraints. 

Such a rapid, out of equilibrium process must result in 
additional entropy production 20 from the latent heat 
of the transition (just like at sudden condensation of su- 
percooled vapour). We start with our model at the line 
where we have quarks with current quark masses, and the 
process ends when we reach the empirical hadron freeze 
out line, where quarks have constituent quark mass. 

It was found by Cleymans et al. [21 that, in a wide 
range of beam energies, detected hadrons freeze out and 
reach the detectors when the average energy of hadrons 
(estimated in a thermal, statistical equilibrium fireball 
model) is between Eh/Nh = (1.0- 1.1) GeV. Just at 
the Freeze Out moment the particles do not interact 
any longer, so the ideal hadron gas mixture is a good 
approximation. Furthermore, if we plot the FO points 
on the temperature-chemical potential [T, ^b] plane, the 
(1.0 — 1.1) GeV constant specific energy contours fall on 
a continuous line indicating that the final state is repre- 
senting a statistical thermal equilibrium state of hadrons 



or a state which is close to it. At the FO state with van- 
ishing interactions and dissipation the expansion of the 
fireball is adiabatic to a good approximation. 

We can also approximate the final FO hadron state by 
assuming a Jiittner distribution to all hadrons. That 
means, we use a relation between chemical potential 
/is and temperature T by means of average energy per 
hadron 



N 

_ h=l 

Nh ^ ^ 

J2 nh{mh,T,fih) 

h=l 



(1.0- 1.1) GeV, (27) 



where hadronic energy density and hadronic particle den- 
sity are given by 



Tcosh 



nh 



9h 
47r2 

grsinhr^U^R'.r!!^^ 



Xr(^ 



mh 
T 



(28) 



Here, ruh is the mass of the hadron, is the degeneracy 
factor, and for baryons we have Hh = I^b while for mesons 
/ift = 0. The Boltzmann- Jiittner gas approach is a good 
approximation to the full Fluid Dynamics and Boltzmann 
equation statistical result: in the medium [T, /is]-range 
the FO contour lines are shown in Diagram (A) of FIG. U 
we have taken into account the lightest 100 hadrons, i.e. 
{h = TfO, 7r±, if", D(1950)). 

In this transition the expansion of the quark gas con- 
tinues until the average energy per hadron reaches about 
(1.0 — 1.1) GeV. In a statistical thermal equilibrium 
model this can be represented on the [T, ^b] plane. In the 
present non-equilibrium model the thermal equilibrium is 
kept, but ^B is not representative of the hadron multi- 
plicity. However, we can approximate the Cleymans-line 
by estimating the hadron multiplicity or density after 
coalescence of our constituent quarks into baryons and 
mesons. The freeze-out curve on the [nB,T] plane is 
plotted in Diagram (B) of FIG. [2) 

The cold uniform nuclear matter has an equilibrium 
density of ub — (1.45 — 1.7) fm~'^, similar to the cen- 
tral density of large nuclei. At higher temperatures a 
considerable number of mesons are also present in statis- 
tical equilibrium, thus the effective baryon density for the 
same energy density decreases with increasing tempera- 
ture. This enables us to construct an effective Cleymans- 
freeze-out line on the [T, n^] -plane for the purposes of our 
non-equilibrium model. 

We have assumed that our rapid hadronization starts 
at the [T, ns] -contour where our quarks have current 
quark mass of about 7MeV. This contour is plot- 
ted in FIG. [3j The chemical equilibrium then ceases 
among quarks and anti-quarks, and we neglect quark 
and anti-quark annihilations in the final expansion stage. 
Thus, the total number of quarks and anti-quarks, Nq 
and N q, remains constant (in contrast to the statisti- 
cal thermal equilibrium model), while the baryon charge. 
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FIG. 2: Diagram (A) shows the freeze-out curves on the 
[^s, r]-plane, corresponding to mean energy Eu/Nu = 
1.0 GeV and 1.1 GeV per hadron. The solid lines were com- 
puted assuming that the hadrons have a Jiittner distribution, 
see Eq. (27 1 (upper solid curve for 1.1 GeV, lower solid curve 
for 1.0 GeV). The dotted lines are from the statistical ther- 
mal model [3] (upper dotted curve for 1.1 GeV, lower dotted 
curve for 1.0 GeV). Diagram (B) shows the same curves trans- 
formed to the [n_B,T] plane. 



Nb — {Nq — Nq) /3, remains also constant as required 
by baryon conservation. 

At the final point of the non-equilibrium expansion 
quarks will coalesce into mesons, baryons and anti- 
baryons, based on the phase space arguments used in 
P5] . The rate of recombination is given by 



— Crr 



gr, 



(29) 



9b 



q + q + q^b : rib = Cb — — — nnUqUq, (30) 

9qgq9q 

_9b 

Hq Tlq Tlq . 

9q 9q 9q 



q + q + q ^ b : Ub = Cb 



(31) 



Here, 



q (q) denotes quarks (anti-quarks), m denotes 
mesons, b (b) denotes baryons (anti-baryons), gm and 
gb are mesonic and baryonic degeneracy factors, and the 
coefficients Cm and Cb should be determined by includ- 
ing phase space factors and an additional normalization 
constant to satisfy the following relations. Notice that 



FIG. 3: The initial conditions for non-equilibrium expansion 
are determined by vanishing quark condensate, see Eqs. Q 
and ([8|. This corresponds to temperature and density con- 
ditions where constituent quark mass equals current quark 
mass Mf — ruf according to Eq. (10 1. 



the same statistical phase-space factors appear in coa- 
lescence rates and in the statistical thermal model if the 
overlap integrals of the coalescing constituents are sim- 
ilar. Eventually other channels via di-quark formation 
can also be taken into account, and the time integrated 
solution will provide the final hadron abundances. 

On the other hand the baryon charge is the excess num- 
ber of baryons over anti-baryons, and this leads to the 
normalization requirements for the final hadrons; 



rib ~ rib ■ 



(32) 



Thus, a fraction "a" of the anti-quarks may form anti- 
baryons: 



Ub 



rib = 



3 ' 3 
and the rest of quarks form mesons 

= (1 - a)nq . 



(33) 



(34) 



The only coefficient, a, should arise from the coalescence 
factors above. Initially from Eqs. (301 - (311 the ratio of 
formed baryons and anti-baryons is rib /rib = Q^, where 
Q = Hq/riq. Assuming that the initial recombination is 
dominant the conservation laws then yield 

Q 



1 



(35) 



The full integration of the rate equations may lead to a 
change of a, which would modify the ratio of anti-baryons 
to mesons to a smaller extent. In [21] the same coa- 
lescence model explained the constituent quark number 
scaling of the flow parameter, V2, from a weak elliptic 
asymmetry of the coalescing quark distributions. 

In this way at any stage of expansion, from an ini- 
tial volume Vb to a point at time t with a given volume 
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V{t), we can get all the meson and baryon densities and 
masses, as well as the total energy from energy conserva- 
tion (neglecting the mechanical work done by the negli- 
gible pressure). 

From the baryon and meson densities ns and n„i and 
the hadronic energy density, E^, we can calculate the 
baryon and meson chemical potentials (the meson density 
may exceed the thermal equilibrium value), and /Xm, 
and the respective temperatures Ti,{t) and Tm{t). Notice 
that with the broken chemical equilibrium the equality 
of chemical potentials is broken. In this rapid hadroniza- 
tion process the thermal, pressure and flow equilibrium 
will also break down. Thus the details how this happens 
should be described in terms of extensive variables as we 
will see in the next Section. 



B. Choice of initial conditions 

Let us consider now the initial conditions for all these 
parameters. The applicability of Eq. (10) is restricted 



and using Eq. [26] to get the chemical potential and sub- 
stituting it into Eqs. p7|, fill), and (flQl). 



by the condition that the constituent quark mass has to 
be positive. Accordingly, we will take the initial condi- 
tions such that in the initial state the in-medium chiral 
condensate {qq)nB,T vanishes and quarks have a current 
quark mass. This condition implies, that the density and 
temperature dependent terms in Eq. ( 10 ) are equal to 
1. The corresponding curve in the temperature-density 
plane is shown in FIG. [3] We used for the current quark 
mass ruf = 7MeV, ignoring the difference between u and 
d quarks. 

In this model it is assumed that the chemical equilib- 
rium between the quarks and anti-quarks breaks precisely 
on this curve, and the quark and anti-quark numbers are 
conserved separately during further expansion. At the 
earlier stage, the quark and anti-quark numbers, Nq and 
Nq, are computed from the temperature T and baryon 



number density ub using Eqs. ([26j), (24) and (25) with 
chemical potentials ^.q = ^b 1^ and JLq — —iib/S. But af- 
ter crossing the line on FIG. [3] the numbers of quarks and 
anti-quarks, Nq and Nq, are assumed to stay constant. 
During further expansion /iq and Jlq change separately 
and can be obtained numerically from the quark and 
anti-quark densities, n„ and n„, by means of Eqs. (24) 



and ( 25 ) . Since this is a non-equilibrium model, all nu- 
merical calculations are carried out on the [nB,T] plane, 
using the densities ub, Uq, riq and the temperature T 
(and not on the [/i^,?"] plane). 



After the initial moment, the quark and anti-quark 
densities decrease inversely proportionally to the vol- 
ume V of the system: rig = Vb/Vn°, Uq — Vo/Vnf^, 
nB — Vo/Vrig, where Vq is the initial volume. Then the 
total entropy of the system can be expressed simply as 
a function of volume, temperature and numerically ob- 
tained chemical potentials (or temperature and baryon 
density) using Eqs. ([T7|, ^ and (l9^ . 

Using the condition that the total entropy is constant, 
i.e. the entropy density also decreases as s(T, V) = 
Sq{T,V) + Sq{T,V) = Vo/V So, the expansion trajecto- 
ries on the [nB,T] plane can be calculated numerically. 
These trajectories are plotted in FIG. |4j 

The trajectories corresponding to constant energy ex- 
pansion (iso-ergic) were calculated, in a similar way, us- 
ing Eqs. (14) and (15). This case corresponds to a dis- 
sipative expansion. It is apparent from FIG. |4] that the 
adiabatic expansion where s = sq ub /n^B l^ads to the 
fastest temperature decrease than the iso-ergic one where 
e = eo nB/ri%. 
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FIG. 4: Expansion of the gas of quarks and anti-quarks. The 
thick solid line is calculated assuming the quark masses to 
be equal to the current quark mass, that means it is iden- 
tical with the curve in FIG. [3] Beyond that curve we take 
Mf = rrij . The dots on the thick solid line indicate the initial 
temperature and initial baryonic density for each curve of adi- 
abatic and iso-ergic expansion. The thin sohd lines represent 
the trajectories of the adiabatic expansion of the gas of quarks 
and anti-quarks. The trajectories of the iso-ergic dissipative 
expansion (dotted lines) are plotted for comparison. 



C. Expansion of the quark gas 

The expansion of the gas of quarks and anti-quarks is 
considered to be adiabatic, i.e. at constant entropy. The 
initial total entropy density of the gas, sqi can be calcu- 
lated from the initial temperature, Tq, and initial baryon 
charge density, n^, assuming chemical equilibrium be- 
tween the quarks and anti-quarks, i.e. /Ug = —JLq = UbI'^, 



D. Recombination into hadrons 

Due to confining forces, the quarks will finally recom- 
bine into hadrons. We assume that this happens rapidly 
at the point of recombination when the average energy 
per hadron (including the background field) decreases to 
Efi/Nu = 1.2 GeV, a value which is still above the val- 
ues of (1.0 — 1.1) GeV obtained by Cleymans et al. [3]. 
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This corresponds nearly to the energy per hadron of the 
empirically observed freeze-out. The endpoints of the 
expansion curves where the recombination happens are 
shown in FIG. [5l 




0.00 0.1 



FIG. 5: A series of curves of adiabatic expansion (thin solid 
lines) of the gas of quarks and anti-quarks. The dots on the 
adiabatic expansion curves indicate where the rapid freeze- 
out and hadronization happens. These points are determined 
based on the condition that the energy of the system, includ- 
ing the background field, divided by the estimated number of 
hadrons reaches 1.2 GeV hadron. The thick solid line is the 
same as in FIG. |3l 

The quarks are assumed to recombine into three types 
of particles: baryons, anti-baryons and mesons, which 
contain three quarks, three anti-quarks, or a quark and 
an anti-quark, respectively. All hadrons are assumed to 
have a mass that is the sum of the masses of their con- 
stituent quarks at the freeze-out, i.e. 



Mh = = 3 M 



FO jiFO 



B 



M^, = 2M 



f 



FO rpFO 



(36) 

(37) 



Further differences between the various hadron species 
are disregarded. Most of the anti-quarks will pair with 
quarks to form mesons, but a small fraction, a, will form 
anti-baryons. This ratio, a, can be estimated based on 
the recombination rates given in |23| . Thus the baryon, 
anti-baryon and meson densities nt, fib and rim are cal- 

and 

331 ' 



culated from the quark and anti-quark densities {uq 



j) using Eqs. 



(jM]) and 

We assume that at the point of recombination the flow 
freezes out too, and both the thermal and chemical equi- 
libria cease. The hadrons are assumed to have Jiittner 
distribution after the freeze out, but the temperature pa- 
rameter in this distribution will be different for baryons 
and mesons. The parameters of the distribution after re- 
combination are calculated from the condition of energy 
conservation: the thermal energy of each hadron type 
will be equal to the energy of their constituent quarks, 
E^^ = El^. Due to the different masses of baryons 
and mesons, their temperature parameters will be dif- 
ferent. The temperature ratio T^/Tm will correspond 
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FIG. 6; Baryon (dotted) and meson (solid) distributions as 
functions of transverse momentum pt (Diagram (A)) and 
transverse energy Et (Diagram (B)), for the case of an adi- 
abatic expansion. The initial state for this calculation was 
0.21 fm"^ and To = 176 MeV. 



to the mass ratio Mb /Mm = 3/2. The distributions 
of baryons and mesons for a calculation done with ini- 
tial state n% = 0.21 fm"^ and To = 176 MeV is shown 
in FIG. [6] and FIG. [7j Here the final baryon and me- 
son temperatures are Tb = 228 MeV and = 152 MeV 
and Tb/Tfn = 3/2, while the final constituent quark mass 
is 308 MeV. The scaled pt and Et distributions become 
identical under this condition, however this is not enough 
to reproduce the NCQ scaling of V2{pt) indicating that 
the recombination influences the flow velocities of the fi- 
nal hadrons. This concept was already pointed out in 
[25] based on the properties of the collision integral. It 
is important to point out that the transport theoretical 
treatment and the collision integral are applicable also at 
situations when the local equilibrium has ceased to exist. 



V. ELLIPTIC FLOW 

The elliptic flow parameter, V2, can be calculated from 
the final, post Freeze-Out distribution by the Cooper- 
Frye formula. Assuming an isochronous FO hypersur- 
face, we obtain simple expressions for final measurables 
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FIG. 7: These diagrams show the same distribution func- 
tions for baryon (dotted) and meson (sohd) as calculated in 
FIG. |6] but re-scaled according to constituent quark number 
ricq {ricq = 3 for baryons, Ucq — 2 for mesons). Diagram 
(A) shows the scaled distribution Ucq ft as function of re- 
scaled transverse momentum ricqPt- Diagram (B) shows the 
re-scaled distribution ricq ft as function of re-scaled transverse 
energy UcqEt- These re-scaled distributions almost coincide 
for baryons and mesons; they are on top of each other so no 
difference can be seen here in the diagrams. The initial state 
for this calculation was Un = 0.21 fm"^ and To = 176 MeV. 



A. Formula for the elliptic flow 

In this Section we mainly follow the arguments of 
Refs. [211 [251 ■ The kinematic average of a quantity 
A{x,p) is given by 



{A) = 



I"'" J 


f d^pf{x,p) A{x,p) 


J 


f J 


f d^pf{x,p) 



(38) 



where f(x,p) is the one-particle distribution function. 
Especially, we are interested in the elliptic flow V2, that 
is the following kinematic average 



V2 



pI-pI 

, Pi 



(39) 



According to Eqs. (38 1 and (39), we obtain 



V2 



d'x I d^pfix,p)?^^ 

Pt 



d^pf{x,p) 



(40) 



Now we subdivide the system into N fireballs or cells, 
each of which has a given volume Vi and contains a given 
number of particles iV^, which are distributed by a given 
distribution function fi{x,p) inside cell i. Since in a 
given cell the elliptic flow parameter V2 does not depend 
on coordinate x, we can take the integral over d^x and 
obtain 



jy 27T OO CO 

^ J d(l) J dp\\ j dpt Pt r{p) cos2(j) 

^—1 n n n 



V2 







N 



l-n OO OO 



(41) 



Ap\\ J dpt Pi f'ip) 







where we have used the relation cos 2 — cos (j) — 

2 2 

sin^ 6 = ^ % with — ^ pi + p?, being the trans- 

Pt Pt y X y ^t, 

verse momentum. In our study we consider the mid- 
rapidity particles, {pz = or the rapidity y = 0), i.e. 
p — (px,py) and there is no integration for the longitu- 
dinal momentum or rapidity. Furthermore, we are in- 
terested in the momentum dependence of elliptic flow 
parameter V2. Thus we obtain from (411 the expression: 



V2 {pt,y = 0) 



N y 

^ d(j) f-ip^^py) cos2(?!) 



N 

^ y d(bfip.,Py) 



, (42) 



where the sum runs over all individual cells i = 1, N, 
and we have used the fact that for mid-rapidity particles 
p^ — Px + pI = Pt- Let us study the elliptic flow of a 
single type of particle. We assume a Jiittner distribution 
in each individual cell of the bag, see Eq. ( 12 ), and obtain 



{2TThy 
1 

(2 7r;i)' 



exp 



exp 



- 7»P0 + 7' VxPx + 7' v'yPy 



V - liPa+Y v'Pv 



(43) 



where p^, is the component of the momentum parallel to 
the cell velocity v*, and /i* is the chemical potential of 
the given particle type in the cell i; here we have also 



used < = 7'(I,<,w*,0) , 
P^ = {Po,Px,Py,Q) , P^i 



{Po, 



, 0) and 



-Px,—Py,0)- Since the 
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elliptic flow parameter is calculated after chemical freeze- 
out has happened, each particle type has their own chem- 
ical potential 



1 



(27r?i)3 """P V T 



ATTMfT,K2{M,/T,) 



(44) 



where rii is the density and Mi is the mass of the given 
particle type. For particles at CM rapidity the ze- 
roth component of four-momentum equals the transverse 
mass, i.e. it is just the energy of one particle 



^Mf+pl = Ml 



(45) 



Then, according to Eqs. M2| and (43), we obtain the 



following expression for the elliptic flow: 

N y 

xiPy) 

-2 = , (46) 



N 



i=l 



X : Py) 



where the sum runs over all cells of the bag. Now let 
us insert Eq. ( 43 ) into expression ( 46 1 . Then, by means 



of Eq. ( 44 1 we obtain the following expression for the V2 



parameter: 

N 



Eiv, 



V2 



i=l 



N 



i=l 

where TV, denotes 



27T 



cos2</) e^'^'P"/^- 



(47) 



T,K2{M,/T,\ 



(48) 



If the direction of the velocity of cell i relative to axis x is 
denoted by 0g, then p„ can be written as p„ = pt cos(0 — 
(/)o), and the following expression is obtained: 



N 



J2 N^ e--^"^^' /T. cog 20j, /2(y V, pt/T,] 



V2{pt 



4=1 



N 



(49) 



Jo and I2 denote the zeroth and second order modified 
Bessel functions of the first kind [26] . 

Notice that if the source temperatures are all equal or if 
the Freeze Out temperature of all fluid cells is the same, 
than one observes a linear increase of V2 as a function 
of Pt or pt/ncq- This is frequently quoted as a linear 
hydrodynamical increase, observed already in the first 
model calculations, e.g. [21] ■ Such hnear increase leads 



trivially to a constituent quark number scaling. It is easy 
to see that this feature is just a consequence of the equal 
temperature assumption. If for example we assume a 
large, hotter, static central source, the high behaviour 
of V2 will be dominated by this hot source and thus V2 
will decrease at high pt as pointed out in ref. [25] . 

If we assume that the temperatures of all cells are 
the same, Ti — T, then the temperature- and mass- 
dependent parts of Ni (see Eq. |48|) cancel from the nu- 
merator and denominator of Eq. ( |49[ ), and wc get (note 
at FO the constituent quark mass does not depend on 
density and temperature anymore) 



N 



V2{Pt) 



Y,N.,e-^^'''^'^ cos2(/)j, hil'v^Pt/T) 

i=l 

N 

E^4e-^'*^'/^/o(y^;,Pt/T) 

i=l 



(50) 



In the special case of four cells moving into the four direc- 
tions (±a;,±y), Eq. (50) reduces to the expression given 
in [24]- 

The simplest configuration that can approximate ellip- 
tic flow is dividing the system into two cells that move 
in opposite directions with the same velocity v. For this 
"two-cell" model case the V2 parameter is expressed as 



V2{Pt 



hjjvpt/T) 
loilvpt/T) 



(51) 



A slightly more complicated possibility is having one 
larger non-moving central cell, and two smaller side cells 
moving in opposite directions, as shown in FIG. |10[ In 
this "three-cell" model case V2{pt) can be expressed as 



V2{Pt) 



2Nse-'^^'''^l2{lvptlT) 
2Ns e-r^'t^/T'Ioi-ivpt/T) + Nc e'^'^^^ 



(52) 



Here denotes the particle number of the central cell, 
while Ns denotes the particle number of the identical 
side-cells. This configuration corresponds to a flow with 
less pronounced asymmetry. The large central cell has a 
spherical momentum distribution, while the smaller side 
cells introduce a slight asymmetry to this. 



B. Calculation of the elliptic flow 

1. Two-cell model 

As a first approximation, the V2 parameter was cal- 
culated using a simple model of elliptic flow where the 
system is divided into two droplets moving in opposite 
directions with the same velocity. For this model, V2 is 
given by Eq. (51). The baryons and mesons were given 
different flow energies, such that the ratio of flow en- 
ergy per quark is {FEh/n^q) /{FE^/ricq) = 3/2. This 



leads to different flow velocities for baryons and mesons. 
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and reproduces the constituent quark number scaling of 
the elliptic flow parameter. In this model the scaling 
according to the transverse momentum ptand transverse 
energy Et are tied to each other, i.e. the scaling is either 
present or not for both these variables. By this point of 
the evolution of our system, all energy in the background 
field is exhausted and the internal, excitation and ran- 
dom kinetic energies of the hadrons reach the FO value 
(1.0 — 1.1) GeV. The V2{pt) curve obtained with the ini- 
tial state parameters discussed in the previous Section is 
shown in FIG. [3 
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FIG. 8: The V2 parameter as a function of pt, calculated 
from the two-cell model, Eq. ( |51[). The initial state used is 
the same as in FIG. [6] and FIGTI?] The dotted curve repre- 
sents the baryons, while the solid curve represents the mesons. 
The cell velocities for baryons and mesons are vt = 0.26 and 
Vm ~ 0.21, corresponding to a flow-energy ratio of 3/2 of the 
constituent quarks of the two different particle types (calcu- 
lated relativistically). 

The elliptic flow parameter, re-scaled according to con- 
stituent quark number is shown in FIG.[9j In this flgure, 
the baryon- and meson-curves coincide for low pt, i-e. 
constituent quark number scaling of V2 is reproduced for 
small values of the transverse momentum, but not for 
Pt > 400 MeV. These results of QNS can be further im- 
proved by considering a three-cell model, as we will see 
in the next Subsection. 



2. Three-cell model 



FIG. 9: The re-scaled elliptic flow parameter, vijucq, as a 
function oipt/ricq (Diagram (A)) and E%/ncq (Diagram (B)), 
calculated from the two-cell model. The dotted curves repre- 
sent the baryons while the solid curves represent the mesons. 
The curves coincide for lowpt value, i.e. the constituent quark 
number scaling is reproduced for the low pt region. 




FIG. 10: Scheme of three-cell model. The asymmetric flow is 
approximated by dividing the system into either two cells that 
move in opposite directions with velocity v, or two smaller 
moving side cells and a stationary central cell. The schematic 
of the latter configuration is shown here. 



The asymmetry in the two-cell model described in the 
previous Subsection is very strong. In [23] it is shown 
that the constituent quark number scaling is more precise 
if the V2 coefficient is small, and the higher harmonic 
coefficients, [k > 2) are negligible. Therefore, a three- 
cell model with one large stationary central cell and two 
moving side cells was also studied. The schematic scheme 



of this arrangement of cells is shown on FIG. 10 



The elliptic flow parameter of the three-cell model is 
given by Eq. ( 52 ) and the result is shown in FIG. 11 The 



side cells were assumed to have the same velocities as in 
the case of the two-cell model, and the temperature of 
all cells was the same. The particle number ratio of the 
central cell to the side cells was set to Nc/Ng = 10. It 
should be noticed, that the elliptic flow parameter, shown 
in FIG.jTT] is not insensitive on such choice, e.g. a ratio 
of Nc/Ns — 5 would enlarge the value of V2. 

The elliptic flow parameter, re-scaled according to con- 
stituent quark number is shown in FIG. |12[ We have 
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(A) 




FIG. 11: The V2 parameter as a function of pt, calculated 
from the three-cell model, Eq. ( |52[) . The initial state used is 
the same as in FIG. [6] and FIG. [7| The dotted curve repre- 
sents the baryons, while the solid curve represents the mesons. 
The cell velocities for baryons and mesons are vt = 0.26 and 
Vm ~ 0.21, corresponding to a flow-energy ratio of 3/2 of the 
constituent quarks of the two different particle types (calcu- 
lated relativistically) . A particle ratio Nc/Na = 10 has been 
assumed. 



found that QNS is insensitive on the chosen particle num- 
ber ratio Nc/Ns = 10, e.g. a ratio Nc/Ng — 5 yields very 
similar results as shown in FIG. [12] According to the re- 
sults in FIG. [12] the three-cell model is able to reproduce 
the constituent quark number scaling of V2 for a wider 
range of pt values compared to the two-cell model. It 
must be noted that this is not only a result of the reduc- 
tion of the momentum distribution asymmetry compared 
to the two-cell model. However, the choice of flow veloc- 
ities for baryons and mesons is still relevant. 



VI. SUMMARY 



(B) 




FIG. 12: The re-scaled elliptic flow parameter, v-ijucq, as a 
function of pt/«-c<j (Diagram (A)) and E%/ncq (Diagram (B)), 
calculated from the three-cell model. The dotted curves repre- 
sent the baryons while the solid curves represent the mesons. 
The curves overlap nearly completely, and the constituent 
quark number scaling is reproduced for a much wider range 
of Pt than in the case of the two-cell model. The velocities of 
the side sources are the same as in FIG. [8] A particle ratio 
Nc/Ns = 10 has been assumed. 



A model of rapid hadronization was considered in 
which the effective constituent quark mass depends on 
density and temperature. In this model the evolution 
starts from an ideal QGP in chemical, thermal and me- 
chanical equilibrium. Then, one after the other the chem- 
ical, thermal and mechanical equilibrium break down 
rapidly, while the quarks build up constituent quark 
mass, and the background gluon field (bag constant) 
breaks up and vanishes. This model can be considered as 
a simple representation of the breaking chiral symmetry 
and deconfinement in a dynamical transition crossing the 
Quarkyonic phase J]. 

The elliptic flow parameter V2 was calculated for the 
final hadron distributions obtained from the model, and 
the constituent quark number scaling was partially re- 
produced. We assume the following stages of hadroniza- 
tion: 1) The chemical equilibrium between the quarks 
and anti-quarks breaks and the chiral symmetry breaks 
at the same time. 2) The quark and anti-quark numbers 



are assumed to be conserved during further expansion, 
i.e. quark, anti-quark pair creations and annihilations 
are ceased. 3) As the quark gas expands and cools, the 
quarks gain mass according to Eq. (10), and the field 
associated with the Bag constant B decreases. 4) The 
thermal freeze-out and recombination into hadrons are 
completed when the mean energy per hadron reaches 
the FO value of Eh/Nh = (1.0 - 1.1) GcV. The cre- 
ated hadrons are not in local thermal and flow equilib- 
rium with each other. We considered only three types 
of hadrons: baryons, anti-baryons and mesons, ignoring 
the differences between different hadron species, and the 
momentum distribution of hadrons was calculated. 

Using the obtained momentum distributions, the ellip- 
tic flow parameter V2 was computed from a simple multi- 
source model of the elliptic flow. Two cases were con- 
sidered. In the first case, the flow asymmetry is approx- 
imated by two fireballs moving in opposite directions. 
Although this leads to a very strong asymmetry in the 
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momentum distribution, it leads to a V2{pt) curve that 
scales with the constituent quark number for small pt 
values only. In the second case, we added a large central 
fireball to reduce the asymmetry and approximate pure 
elliptic flow better (i.e. reduce the higher Vk {k > 2) 
harmonic coefficients) . This three-source model was able 
to reproduce the constituent quark number scaling for a 
wide range of pt values. 

The present model is highly simplified and attempts to 
provide an insight to the rapid hadronization and freeze 
out process in view of the quark number scaling. We con- 
sidered the features arising from breaking down of equi- 
librium in this process in terms of thermo and fluid dy- 
namical parameters which are applicable to partial com- 
ponents of the matter. We intend to implement these 
concepts in more complex models (like hybrid models) 
and to search for more fundamental reasons for the ob- 
served freeze out features in terms of the partial exten- 
sives. 

Despite the simplicity of the model, it is capable of 
reproducing the constituent quark number scaling of the 
V2 elliptic flow parameter. The presence of constituent 
quark number scaling in experimental data suggests that 
the elliptic flow develops in the Quark-Gluon Plasma 



phase, before the quarks recombine into hadrons. There- 
fore understanding the origin the elliptic flow can provide 
insight into the quark phase of matter. 
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